Entangling capacities of noisy two-qubit Hamiltonians 
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We show that intrinsic fluctuations in system control parameters impose limits on the ability of 
two-qubit (exchange) Hamiltonians to generate entanglement starting from mixed initial states. We 
find three classes for Gaussian and Laplacian fluctuations. For the Ising and XYZ models there are 
qualitatively distinct sharp entanglement-generation transitions, while the class of Heisenberg, XY, 
and XXZ Hamiltonians is capable of generating entanglement for any finite noise level. Our findings 
imply that exchange Hamiltonians are surprisingly robust in their ability to generate entanglement 
in the presence of noise, thus potentially reducing the need for quantum error correction. 
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Considerable experimental efforts have been devoted 
in the past few years to the creation of entangled states, 
with impressive success in systems such as trapped ions, 
coupled atomic gas samples, polarized photons, and most 
recently, superconducting qubits lj. An important mo- 
tivation comes from quantum information processing 
(QIP), where entanglement is believed to play an impor- 
tant role in algorithmic speedup, communication tasks, 
and cryptographic applications [2|. As the generation 
of entanglement often involves the manipulation of an 
interaction Hamiltonian, recent theoretical work has fo- 
cused on the entanglement capabilities of such Hamiltoni- 
ans. In particular, questions concerning optimality 
equivalence classes ||, and entangling power/capacity 
[g, have been raised and answered, under the assump- 
tion of noiseless controls. Here we take the first step 
toward addressing what happens when this assumption 
is relaxed. In particular, we wish to find out the answer 
to the following question: What are the limits imposed 
on entanglement generation via two-body Hamiltonians 
by fluctuations in system-control parameters? [?]. We 
note that, as is well known, quantum error correction 
offers a solution to both decoherence and the type 
of control errors we consider here; however, this solution 
involves a high cost in extra qubits and logic gates. In 
view of the central importance of entanglement in QIP, 
it is of significant interest to find out the limits imposed 
on entanglement generation via interaction Hamiltonians 
and state preparation, without any error correction. 

The Model. — Almost all quantum computing propos- 
als are governed by interaction Hamiltonians that are 
used to enact two qubit operations. The most general 
two-qubit ( "exchange" ) Hamiltonian has the form 



H 
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where S a 



\a a are the angular momentum operators 



(a a are the Pauli matrices), and i, j are qubit indices. By 
applying local unitary operations it is always possible to 
transform H into a canonical, diagonal form |3(. Hence 
we will restrict our attention to the case J^g = J^S a p 



from now on. The various models are then classified as 



follows: XYZ: J x ^ J y ^ J z , XXZ: J x = J, 



= 0, Heisenberg: 



Jx Jy 



+ Jz, 

, Ising: 



XY. J X Jy,J; 

Jx — Jy — 0- 

Two qubits can be entangled by first preparing a prod- 
uct state and then running the interaction for a de- 
sired amount of time, to generate, e.g., a C-NOT or C- 
PHASE gate 0. Tunability of the coupling constants 
in H need not always be possible, even though it is a 
common assumption of QIP proposals. This leads to two 
qualitatively distinct scenarios we consider in this work: 
(a) Tunable interactions - where the interaction can be 
switched on and off (e.g., an exchange interaction medi- 
ated by a tunable tunneling barrier (b) Non-tunable 
interactions - where the interaction is always on (e.g., 
a Coulomb interaction 0), thus requiring, e.g., external 
single qubit operations to refocus the interactions and 
enable controlled entanglement generation. Recent work 
has addressed the problem of universal quantum com- 
putation with non-tunable couplings [l0|. and even un- 
known parameters In the laboratory, however, the 
execution of every single and two qubit operation will 
generally be noisy due to system and experimental imper- 
fections, over which we have limited control. We consider 
phenomenological noise models to describe noisy single 
qubit and two qubit operations, wherein certain control 
parameters vary stochastically. Specifically, we have con- 
sidered two models: (a) Gaussian and (b) Laplacian pa- 
rameter fluctuations. The Gaussian model has universal 
applicability in the case of noise due to many weakly 
coupled random sources (by the central limit theorem). 
It has been extensively used and discussed in stochastic 
quantum mechanics (e.g., [l^L We consider the Lapla- 
cian model mainly to test the robustness of our results. 
Another important model is 1// noise due to bistable 
random fluctuators, which will be considered in a future 
publication. 

Noise model. — Given a Hamiltonian K(t, J) (where J 
is a parameter or set of parameters), a unitary transfor- 
mation U is generated by evolving under K for some time 
t: U{4>) = Texp(— i JJ K(t, J)dt), where T denotes time 
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ordering. In our analysis below we only deal with piece- 
wise constant Hamiltonians; then the angle(s) cf> = tJ . 
An initial state p transforms as p ~ * p(<fi) — U (4>) pU {<f>y . 
We now assume that <f> is Gaussian distributed with mean 
4> (the desired angle) and standard deviation (s.d.) A: 
<f> ~ N(cf>,X). This may be the outcome of Gaussian 
noise in J, r, or both. Thus, under a noisy control the 
actual transformation is 

-\ 1 f°° (*-f) 2 
P -> Pnoisy (A, 0) = ~j=^ J e p((f>)d(j), (2) 

Below we take K to be either an exchange Hamiltonian 
with noisy coupling constants J l J (in which case we as- 
sume for simplicity an equal s.d., denoted CI), or a noisy 
single-qubit Hamiltonian needed to refocus an always-on 
exchange Hamiltonian (and denote the s.d. by A). 

Case (i): Tunable Ising interaction. — In the Ising 
model, the exchange Hamiltonian takes the form Uzz = 
\Ja\a\ = JSlS%. Consider the preparation of a max- 
imally entangled state, starting from the initial state 
Pi = 1 00) (00 1, in the presence of noisy interactions 
|13| . Without noise, application of the Hadamard 
transform Uh(jt) = exp(— iirSy) on both qubits, fol- 
lowed by U Z z(k) = exp(-i J Q T \ Jo\o\ db) = e" 1 ^^ , 
prepares the maximally entangled pure state |£) — 
\ (| 00) + i |01) + i 1 10) + 1 1 1) ) . In the presence of noise, 
the action of both Uh and Uzz must be averaged over a 
distribution of angles, as in Eq. J2Jl . The noisy Hadamard 
transform is a rotation about the y-axis with average an- 
gle 7r and s.d. A, resulting in the mixed state po (A) = 



/2tt\ 



(U H (9)piU H (e^) d9, which can be easily 



evaluated. The mixedness of po is measured by its von 
Neumann entropy: M (p) = — Tr[plog 2 p], as a function 
of the noise parameter A. M (p) — for a pure state; 
M (p) = 2 for a maximally mixed state. In the present 
case the entropy rises rapidly from zero (at A = 0) and 
reaches its maximum of 2 for A ~ 2. Next we apply the 
noisy version of the Uzz gate, with angle Jt ~ N{%, CI). 
The resulting density matrix p (A, CI) is easily computed 
but is not particularly illuminating; instead we present 
the result of using the partial transposition test for entan- 
glement 01 : a 2 (g) 2 state is entangled iff it has negative 
partial transpose (NPT). One then arrives at the follow- 
ing condition for inseparability for the density matrix 
p(X,Cl) 



\2 1 / \ 2 I CI 2 \ 



(3) 



Figure ^ illustrates this condition. Observe that there is 
a significant region of entanglement in parameter space, 
with a trade-off between the tolerated level of noise in 
state preparation and interaction. Interestingly, except 
along the cut A = 0, the transition from entangled to 
separable is sharp. Solving the inequality J3J, we find 
that the condition for entanglement is 




A<{-21og[(e- 



1) 



1/2 



1} 



1/2 



(4) 
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FIG. 1: Condition for entanglement as a function of interac- 
tion error £1 and preparation error A in the Ising model. Plot- 
ted is the PPT criterion 10. The transition line between en- 
tanglement and separability is clearly visible, with the region 
above the horizontal plane correponding to entanglement. 



The finite range of the state preparation parameter A 
indicates that the purity of the initial state is crucial. 
However, with a high quality interaction a significantly 
mixed initial state can be tolerated: From the above fol- 
lows that if the interaction is perfect (CI = 0) then A max 
= (-2 In (\/2 - l)) 1/2 = 1.327, meaning that even an 
initial mixed state with entropy 96% as high as the max- 
imally mixed state would still enable entanglement gen- 
eration. Conversely, if the initial state preparation does 
not involve any noise (A = 0), then the interaction, no 
matter how noisy will be able to produce some entangle- 
ment. Finally, note that for 2 £g) 2 systems of the type we 
are considering here, if a state is entangled then it is dis- 
tillable as well, i.e., one can extract pure states from such 
noisy states using local operations and classical commu- 
nication |l5l| . In particular, Eq. Q therefore guarantees 
that a state is useful for teleportation and all other QIP 
primitives. 

Case (ii): Untunable Ising interaction. — Now we 
do not assume the ability to switch the interaction 
off (as, e.g., in NMR). It is therefore necessary to 
refocus the interaction using single qubit operations 
|l6|. This is done by pulsing an external magnetic 
field along the x -axis (we choose the first qubit for 
this operation). We assume that such pulses can be 
made very fast and strong compared to the interac- 
tion ^(J . Formally, let A, Y, Z be operators satisfying 
su(2) commutation relations: [A, Y] = iZ and cyclic 
permutations (e.g., the angular momentum operators 
S a ). Then it follows from the Baker-Campbell-Hausdorff 
formula that upon "conjugation by ip" : o X = 

exp(— iipZ) X exp(iipZ) = Acosy-1- Fsin^. Note further 
that Ue A U^ — e UAu1 for unitary U and arbitrary A. 
Thus: C^oexp(i8X) = exp(-itpZ) exp(i9X) exp(itpZ) = 
exp(i0(A cos if + Ysinip)). One application is "time- 
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FIG. 2: Plot of the inseparability condition @. Entangle- 
ment in the noisy tunable XYZ model is the region where |z| 
is above the surface shown and < 1. 

reversal" , which results from conjugation by it: 

d o e MX = e-^-e^e 1 ^ = e' l9X . (5) 

We then have p(r 2 ) = U r(tt) pqU , where Ur(tt) = 
[C% o exp(-z/; 2 JSlSldt^eM-ifo 1 JS\^dt) = 
e~ % ~ a z a * , which in conjunction with our Hadamard state 
preparation yields the desired C-PHASE gate. The 
parameters must satisfy the condition J{2t\ — r 2 ) = 
7r. In the noisy scenario the recoupling step is imple- 
mented with a rotation around the x axis by an an- 
gle 6 ~ N(n,A) (where A 2 = A 2 + A 2 ., with A, the 
s.d.'s of the independent Gaussian random variables JTj, 
i = 1,2, J fixed); the only change is then Ur(tt) h- > 
Ur(9), and the final state is given by: p(j2,X, A) = 

f e ~^~ (U R ,(9)p U R (0y) d9, which can be an- 
alytically evaluated. The inseparability condition, ob- 
tained using the partial transposition criterion, yields the 
same form as the tunable case [Eq. Q], provided we re- 
place n/2 by A. Note that this is by no means an a pri- 
ori obvious substitution: Q is the two-body interaction 
strength error, whereas A is the error in the single-body 
coupling parameter. 

Case (Hi): Tunable XYZ Hamiltonians. — We now 
consider the XYZ model H XYZ = £ Q =x,y,* J a S^, of 
which, in the noiseless case, the XXZ, XY and Hciscn- 
berg models are special cases. As in the Ising case we 
assume the initial state is pi = 1 00) (00 1 |lj|. In the 
noiseless scenario, we first apply Ux(k) = cxp(—iTrS x ) 
on the second qubit, yielding po — 1 01) (01 1 . This is fol- 
lowed by UxYz(0x,0 v ,(t)) = exp(-i f Q Hxyzdt), where 
9 x . y = Jx,yT, </> = J 2 t. Letting 6 X + 8 y = j and leaving 
<p arbitrary, this prepares the maximally entangled pure 



state |£) = (i 1 01) + 1 10) ) . Thus in the noiseless case, 
and for all the exchange models considered here, there 
is no dependence on J z . In the noisy scenario, we first 
apply Ux(oj) — exp(—iuiS x ) where w ~ N(ir 7 A). Then 
po (A) is a mixture of the states |00),|01) with respective 
weights [1 =F exp(— X r)]/2. Note that now, in the worst 
case scenario (A — > 00) po can only be 50% as mixed 
as the totally mixed state, whereas in the Ising model 
above po could be fully mixed (or 96% mixed if we insist 
on entanglement generation). 

The unitarily transformed density matrix is given 
by po(A) h-> Uxyz {9 x ,9 y , tf>) po (A) U^ YZ (d x ,6 y ,(j>) = 
p{\ 1 9 + ,6~), where ± = X ± 9 y , whose explicit form 
can be found without much difficulty, and again is inde- 
pendent of <f> (i.e., J z ), now even when the initial state 
is noisy. We next integrate over 9 + ~ iV(7r/4, f2) and 
9- ~ N(fF,Q), where Q = (O 2 + fl 2 ) 1/2 , and where O; 
are the s.d.'s of 9 L . Upon applying the partial transposi- 
tion criterion we find the condition for inseparability 

\z\>[l/b*-(l + af/(l-a) 2 } 1/2 , (6) 

where a = e~^ x ~ , b = e~ 2n2 , z = cos(26*~). This condi- 
tion, plotted in Fig. [2J depends on the noise parameters 
X,fl (as in the Ising case), but also periodically on the 
(mean) coupling constants through the distance between 
exchange models 0~ — (J x — J v )t. From < \z\,a, b < 
1, we have the following sufficient conditions for insepa- 
rability [RHS of Eq. © < 0]: b > (l-a)/(l+a) (curve at 
\z\ = in Fig. [ a, and sep arability [RHS of Eq. © > 1]: 
b < (1 - a)/^2(l + a 2 ) (curve at \z\ = 1 in Fig. EJ). 
Thus we find that that the XYZ model with J x ^= J y , 
similarly to the Ising model, exhibits a sharp entangle- 
ment/separability transition as a function of preparation 
and interaction noise. However, in contrast to the Ising 
model, in general there is also a dependence on a third 
parameter, the (noisy) interaction distance J x — J y , so 
the XYZ model belongs to a distinct class. 

Case (iv): Tunable XY, XXZ and Heisenberg 
Hamiltonians. — Note that because of the integration 
over the distribution of J x — J y above we cannot specialize 
to the XY, XXZ, Heisenberg models (J x = J y =; J case). 
Repeating the calculations above now with J ~ N(j, f2) 
we find that the final state p(X, O) is a mixture of the 
states |00) and |* ± ) = ^= (|01) ± |10>), with respective 
weights i(l- e - A2 /2) and i( 1 + e -iA 2 )( 1=Fe -o 2 /2). Such 
states are always entangled as long as the proportions of 
the states |vf ,± ) are different, which clearly holds in our 
case for A, 51 < 00. Hence the final state in the tunable 
XY, XXZ, and Heisenberg models is always entangled for 
all practical purposes. This then is a third class. 

Case (v): Untunable XYZ, XY, XXZ and Heisenberg 
Hamiltonians. — Since the interaction is always on, we 
need to apply external single qubit refocusing operations. 
This is done by pulsing an external magnetic field along 
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the z-axis (we choose the first qubit for this operation). 
After a lengthy calculation we find from the partial trans- 
position criterion, for the XYZ model, the inseparability 
condition 

(A + Bp) 2 (1 - nf (C M - Df (1 + vf < (7) 

where p — e~ A I 2 ,r\ = e~5 A J A — cos (3 cos 8, B — 
sin (3 sin 6, C = sin (3 cos 5,D = cos/3sin<5, /3 = Aw/2,5 = 
Att/4, and A = 4 : Vr i - The XY, XXZ and Heisenberg 

<->x~i>'y 

models ( J x = J y ) corresponds to B = C = D = 0. Inter- 
estingly not only this is achieved when J x = J y but XYZ 
mimicks XY, XXZ and Heisenberg also when A = An, n 
an integer. It can be shown that the inequality Q is sat- 
isfied as long as A, A < oo. Hence in all these exchange 
models the final state is always entangled for all practical 
purposes. In contrast to the tunable case, the untunable 
XYZ, XY, XXZ and Heisenberg Hamiltonians all lie in 
the same class. 

Laplacian fluctuations. — As a test of the robustness of 
our Gaussian-model based conclusions, we have repeated 
the above analysis when the flucutuations in control pa- 
rameters obey a Laplace distribution: under a noisy con- 
trol the actual state transformation is now 

-% 1 f°° 

P -> Pnoisy (W, 4>) = — J e 2" p ((f)) d(j). (8) 

While we find quantitative differences compared to the 
Gaussian model, the qualitative behavior is identical, 
thus bolstering the universality of the classes obtained. 
To illustrate this fact let us reconsider the Ising case. Let 
u> = A (Jl) be the s.d. of the noisy parameter controlling 
state preparation (interaction) in both the tunable non- 
tunable case. We find, repeating the procedure above, 
that the condition of inseparability can be given for both 
the tunable and non-tunable cases as 

4A 2 (ft 2 + 2A 2 + 2X 2 n 2 ) < 1. (9) 

Recall that, similarly, in the Gaussian case the insepara- 
bility condition of the tunable case can be obtained by ap- 
propriate rescaling of the s.d. Further, note that as in the 
Gaussian model, when the initial state is perfect (A = 0) 
the state remains inseparable for all practical purposes. 
By solving inequality © the inseparability condition can 
be written as A < §[((l+(l + ti 2 ) 2 ) 1 ' 2 -^ 2 )/ (l + Q 2 )}^ , 
which shows that, again as noticed in the Gaussian case, 
the purity of the intial state is crucial and actually even 
more so in the Laplace case. Indeed, if the interac- 
tion is perfect , we obtain the inseparability condition 
A < \A/8 = 0.5946, in contrast to the Gaussian case 
where the the threshold value of A in the no-noise inter- 
action scenario is as high as 1.327. 

Conclusions. — The sharp transition found in the Ising 
and XYZ models is reminiscent of the thermal entan- 



glement transition [18j, and suggests an interesting av- 
enue for further research. The surprising robustness of 
entanglement to noise bodes well for quantum informa- 
tion processing with reduced demands on error correc- 
tion. Another interesting implication of our work con- 
cerns entanglement verification: Knowing the underly- 
ing two-body interaction Hamiltonian and corresponding 
level of control, an experimentalist can confidently char- 
acterize the degree of entanglement his/her system can 
generate, without needing to perform a direct, and often 
difficult, measurement of entanglement 0. 
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